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Recent developments in the three-dimensional aerodynamic theory of 
inviscid flow in transonic axial compressors are reviewed briefly. Empha- 
sis is placed on the newly completed lifting surface theory of a transonic 
ducted rotor. The relationship between the lifting surface theory and 
axisymmetric through-flow theories of turbomachines is illustrated; a few 
examples of the additional information obtainable from the new theory 
are then given. Quasi-two-dimensional cascade theory can also be ex- 
tracted from the present analysis and the relevance of cascade theory to 
the actual three-dimensional problem assessed. Details are reported else- 
where, but some of the qualitative conclusions are discussed here. Even 
moderate departure from uniform spanwise loading of the rotor blades, 
for example, leads to a rather profound influence of the downstream 
wakes, suggesting the need for considerable care in applying cascade 
data on a direct quasi-two-dimensional basis. 


The inviscid, three-dimensional, compressible flow through an axial 
compressor rotor or ducted fan can be described in terms of the perturba- 
tion of the incoming flow by the rotor and its wake. If the incoming flow 
is sufficiently uniform and regular, and if stator interference can be 
neglected as a first approximation, then the flow is steady in coordinates 
fixed in the rotor. If, moreover, the perturbations induced by the rotor are 
“small” they can be described by a velocity potential which satisfies the 
convected wave equation. 


1 A major part of this work was supported by the United Aircraft Corporation Re- 
search Laboratories while Dr. Okurounmu was a member of the Research Staff at 
that Laboratory. A portion of the work of J. E. McCune was also supported by the 
Pratt and Whitney Division of United Aircraft Corporation. 
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156 THEORETICAL PREDICTION OF FLOWS IN TURBOMACHINERY 

For steady flow in rotor coordinates, when the relative Mach number 
is everywhere small, the governing equation for the velocity potential 
reduces to Laplace’s equation. In this paper, however, we will be con- 
cerned with transonic rotors (i.e. , rotors operating with subsonic axial 
Mach numbers), while the relative Mach numbers at the tip may be 
supersonic. The relative Mach number at the hub is usually subsonic and, 
therefore, the flow will generally be of a mixed type. The governing 
(linear) equation changes from elliptic to hyperbolic type at the “sonic 
cylinder,” r = r„, where oj 2 r, 2 + U 2 = a 2 . Because of three-dimensional effects, 
however, the linear theory does not exhibit the degeneracy for relative 
Mach numbers approaching unity which occurs in two-dimensional 
theory. 

The coordinate system we use in this paper is fixed in the rotor (fig. 1) ; 
G 3 is the angular velocity of the rotor, U the (purely axial) velocity far 
upstream, a the (undisturbed) speed of sound and r the radius; x is the 
axial coordinate, 6 the azimuthal. The corresponding dimensionless 
variables are ( z,6,a ) = (ux/U ,d,o>r/U) . 

In this paper we will be concerned primarily with the lifting problem. 
The thickness problem was treated earlier (see refs. 1 and 2). In those 
papers, source-type singular solutions ( B radial source “spikes”, where 
B is the number of blades in the rotor) were constructed by superposition 
of the acoustic eigenmodes of the system in a straight annular duct of 
infinite extent. These source spikes, of arbitrary strength, Q(r), were then 



Figure 1. — Geometry and coordinate system fixed in rotor, a is the angular speed of the 

rotor relative to the shroud. 
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distributed over helical surfaces of finite chord approximating the blades, 
thus representing blades of arbitrary thickness distribution (with the 
limitation, of course, that the thickness-to-chord ratio remain small). 
The camber line of the blades was left free in that study so that the 
“rotor” produced no lift and hence no positive work. Attention was thus 
focused on the effects of blade thickness in three dimensions. 

The conditions under which quasi-two-dimensional cascade theory 
could be recovered from the three-dimensional theory were delineated in 
reference 1, while it was shown in reference 2 that large departures from 
quasi-two-dimensionality are to be expected in the transonic regime 
because of an acoustic resonance of certain eigenmodes. It was also shown 
that the wave drag due to thickness should be considerably smaller than 
that corresponding to quasi- two-dimensional (strip) theory. 

Reference 3 describes experiments on a “free-wheeling” transonic rotor 
undertaken in part to verify the latter theoretical result. In this regard, 
the experiment was somewhat inconclusive, partly due to the difficulty of 
identifying and separating out the various types of drag and/or losses and 
partly due to the fact that any actual (rigid) rotor, of course, produces 
lift, varying from hub to tip. For a free-wheeling rotor (zero torque 
input) the result is that vorticity is shed downstream, adding substan- 
tially to the drag. 

In order to complete the three-dimensional potential theory, the basis 
elements of the theoretical lifting problem were set up in reference 3. 
Following the general procedure of superposing the appropriate eigen- 
modes to construct singular (Green’s function) solutions, B bound vortex 
spikes were constructed, having arbitrary strength r(r) . (The method is 
almost identical to that used in ref. 4, except that a finite hub-to-tip ratio 
was included and an important error occurring in ref. 4 was removed.) 

The most important new feature of the lifting problem, relative to the 
thickness problem, is the necessity of including the downstream wakes 
of shed vorticity (one helical wake for each blade) with strength propor- 
tional to dr/dr, the rate of change of bound circulation along the span. 
This is done by a slight modification of the method of Reissner (ref. 5), 
to allow for the presence of the hub and shroud. Thus, a “wake potential” 
is included in the downstream flow, added to the acoustic eigenmodes, and 
construction of the bound vortices at the blades proceeds as before, with 
the wakes now included. It is interesting to observe that the wakes them- 
selves excite acoustic modes, except when r(r) = constant. 

The acoustic eigenmodes mentioned above are versions of the familiar 
“spinning” modes associated with cylindrical geometry (refs. 6 and 7). 
If the relative Mach number at the tip exceeds unity, some of these modes 
propagate undiminished in strength (in the inviscid, linear theory) 
upstream and downstream, while the remaining modes die out exponen- 
tially at large distances from the rotor. The propagating modes are said 
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to be above “cut-off”; modes near cut-off can set up the acoustic reso- 
nances mentioned above. One effect of such resonances is to create a 
significant span wise flow, yielding strong communication between hub and 
tip. Under these conditions (which require supersonic tip Mach numbers) , 
a hub section with a subsonic relative Mach number can nevertheless 
have a pressure distribution over the blade which is more typical of 
supersonic flow than subsonic flow (refs. 1, 2, and 8) . 

In reference 9, the lifting theory was refined and extended. In that 
paper, the blades were characterized by their total bound vorticity, 
r(r), and the first-order static pressure rise across the rotor, the turning 
angles, the torque required, etc., were determined as functions of r. For 
example, the first-order static pressure rise was found to be 


<p.)-<p_.) 2 or / 2 f 1 r(„)\ -2<r r T 

~ 0 * \l-h*J h V V UL T ) p UL T 


( 1 ) 


where 7 j = r/r r = aAr, h = r H /r T , < > denotes the azimuthal average of a 
given quantity, / 3 2 = 1 — U i /a 2 =l—M 2 , and Lt = 2^tt/B, the blade 
spacing. The subscripts T and H denote tip and hub, respectively. The 
sign convention on r is such that it is negative for lifting blades. The 
corresponding torque required, to first order, is 

T<d= — p _ x U 2 r T 2 B [ vdrjT (2) 

' h 


By extending these calculations consistently to second order (a procedure 
analogous to the computation of induced drag and wave drag in ordinary 
wing theory) , the losses due to energy stored in the wakes of shed vorticity 
(when T(r) ^constant) as well as those due to acoustic radiation (for 
supersonic relative tip Mach numbers) were estimated. These are pre- 
sented in figure 2 for a typical rotor (B = 40, h = 0.8), in terms of a 
dimensionless efficiency decrement. Since these results are essentially 
“integral relationships” (i.e. , obtained from momentum balances, etc.) 
we expect them to carry over without change to the lifting surface theory, 
which will be the main subject of our discussion. Despite the linearizing 
assumptions inherent to the theory, we estimate that these results will be 
accurate up to static pressure ratios across the (single) rotor of about 1.3 
(see eq. (1)). It should be noted that the “concentrated bound vortex 
solution” of reference 9 is not a lifting-line theory in the sense of Prandtl, 
since no quasi-two-dimensional assumptions were used. 

The lifting surface theory, in analogy to the procedures described in 
reference 1, can be constructed from the concentrated bound vortex 
solution by distributing the bound vorticity, with its associated wakes, 
over helical surfaces representing the blades. Details will be made avail- 
able shortly in reference 10. In the next sections, we describe the salient 
features of the theory and some of its more interesting results. 
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6= 40, h = 0.8, C 0 „/L T = 1.06 
uT7= a+bi7,D = <r T -r H )/r 
FOR M = 0.6 Os = 1.333 


M = 0.5 os =1.732 



Figure 2. — Losses for a typical 
rotor configuration as caused by 
( 1 ) energy stored in wakes 
(D 7^0) and (2) acoustic radia- 
tion (<rT>tr a ). Solid lines ( left- 
hand scale) give values of the 
efficiency decrement, normal- 
ized to Cl, an average ( over the 
span) sectional lift coefficient. 
The relation between Cl and 
f , defined in the text, is shown, 
for this rotor geometry, by the 
broken line ( right-hand scale). 
The sharp increase in losses due 
to acoustic radiation as <tt ex- 
ceeds a, is readily observable 
for both D =0 and D =0.221S. 
In the latter case, the efficiency 
decrement associated with non- 
uniform loading is comparable 
to the acoustic radiation loss. 
(Figure taken from ref. 9.) 


LIFTING SURFACE THEORY 

The velocity potential obtained by distributing bound vorticity, in the 
manner described, over the helical surfaces representing the rotor blades 
can be written 

where £ labels the axial location of the bound vortex filaments while 


H 

III 

0<Z<C ax (max) 

o 

III 

x<0 

£ = Cax (max) 

X>c ax (max> 


If we denote by <p(r) = tan -1 a = tan -1 (wr/U) the complement of the 
stagger angle at each radius, the axial projection of the local chord is 
Cax(r) =c cos <p, and in (3) c 8X (max) is the maximum of this quantity along 
the span. The elemental potentials bg>' 1 and 5 <t> d are, respectively, the 
upstream and downstream velocity potentials associated with B radial 
strips of concentrated bound vorticity located on the helical surfaces 
x = £, d = w£/U-\-hr/B, Z=± 1 , ±3, ±5, .... Each strip has elementary 
chord d£/ cosy> and strength y (ri.fi) d!)/ cos g>. Their associated wakes of 
free vorticity, proportional to d/dt](y/cos <p), are included in 54> d . The 


Cax (m “ X> \ 

'■da+ ( ' srdz) 

J X ' 


(3) 
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total bound vorticity (per blade) at each radial station is the local blade 
circulation : 


r(r) = T(jjrr) = / 7 (»J,£) d£/cos <p 

J n 


(4) 


The elemental potentials occurring in (3) can be expressed as follows 
(ref. 10; real part implied) 
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( 6 ) 


In these expressions, the “radial eigenfunctions” K b b (*»*»■/ or) —RnB(nnkTi) 
are orthonormal combinations of Bessel and Neumann functions whose 
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properties, including the radial eigenvalues K„k, are determined by the 
requirement of vanishing radial velocity at hub and shroud. Details are 
available in references 1 and 3. The quantities \ nk , determining the axial 
behavior of the acoustic modes, are given by 


nB j K„k~ orM- 


( 7 ) 


Xot = 


KO k 

far 


These quantities are either pure real (yielding axial decay) or pure 
imaginary (yielding propagating, “spinning” modes). Since the K nk are 
all greater than | nB | (ref. 11), the latter possibility occurs only if 
<tt>M/( 3 = ccr s / U, and then only for a finite number of modes at each n; 
k= 1, 2, . . . k„*. The condition a T > a, implies supersonic relative Mach 
numbers at the tip. 

In (6), the variable f =d-ux/XJ = 6-z is the helical coordinate used by 
Reissner (ref. 5) and the (generalized) function is a “sawtooth func- 
tion,” the essential properties of which are that d£i/dd = 1, d£i/dz=— 1 
everywhere, while the function itself is discontinuous (by an amount 
^ir/B) at f = ±ir/jB, ±?>ir/B, .... The combination of the second and 
third terms in (6) makes up the wake potential mentioned earlier, 
representing the free vorticity shed in the wakes and the induced flow 
between them. The wake functions, Xniv,!;), can be written in terms of 
modified Bessel functions of the first and second kind: 


Xn(»J,£) =Xnn{nBr\OT,k) =a„(|)/„B(nS<r) -\-(in{£)KnB{nB<T) 


+J»s(nB<0 j — K nB (nBa')a' — 


d \ y{<r'/<r T ,V) 
cos ip(a') . 


da' 


— K nB (nBa) f — I nB (nBa')a' — 
da da 


d y {a'/<TT,£) 

_ COS <p(a r ) . 


da' 


( 8 ) 


where a„(f) and /3 n (tj) involve definite (radial) integrals over 
d/da[y(y,g)/cos <p~\ and depend on the parameters a H and a T , as well as 
nB (refs. 3, 5, 9, and 10). Note that each X n(n,£) vanishes when t/cos <p 
is independent of o-(or rj) ; i.e. , when T= constant. 

The quantities -y (£) , y nk(£) and h„ k (^) are the coefficients of expansions 
of the functions 7 ( 17 , £) /cos v and x» (>?,£) in terms of the orthonormal 
radial functions R„b («»*!?)■ Thus, 
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y(v,t) 

COS <p(ti<tt) 

y(v£) 


= y(£) + '52yok(%)Ro(nokn) 

k- 1 
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= ^ ) v«A: (£) RnB ( Knkw) , Ti^ 0 
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00 

Xn(v,£) = ^h n k (£) RnB ( Knkv) 


( 9 ) 


The usual Fourier-Bessel formulas for the expansion coefficients apply; 
for example 
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(£) = f v di) 

J h 
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Note that 
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y(v,$) 

COS <p(t]<Tt) 
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(max) 
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as defined in equation ( 1 ), and that the coefficients y 0 *(£), as well as the 
Kk(£) vanish whenever 7 /cos <p is independent of radius. 

In the special case for which c 8X = constant and 7 (t 7 ,f)/cos <p is fac- 
torizable (i.e., 7 /cos <p= T(n)gH ) ) we have the especially simple relation- 
ships, y»*($) = h nk (£) =H nk g($), 7 (() = rg(Z), where, from (4), 

/ c ax 

<K?)#= 1 (11) 

T„k, H„k and T are then identical to the corresponding quantities occurring 
in the concentrated bound vortex solution (refs. 9 and 10) , in which the 
blades are characterized by T(r). 

The three-dimensional lifting surface solution, described formally by 
equations (3) through (10), has the desired property of producing dis- 
continuities (at the helical surfaces representing the blades) of the velocity 
component parallel to the blade surface. In fact, if we denote this com- 
ponent by tv, the discontinuity in tv is — y(v,x), and the blade loading is 
— p_ aj C/ r 7 (r;,x), where U r = (t/ 2 +o; 2 r 2 ) 1/2 . The detailed proof that (3) 
with (5) and ( 6 ) has this property is available in reference 10. At the 
same time ty> the velocity component perpendicular to the helix of 
advance, is continuous, and v V '/U r defines the slope of the blade camber 
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line (in the absence of thickness effects) for any loading. As is usual in 
wing theory, if the blade shape is specified and the loading required, 
the mathematical problem becomes one of inverting an integral equation — 
in our case an unusually complicated one. 

Actually, if one is given a shape for the rotor blades (desired thickness 
and camber) , the thickness problem must be solved first (refs. 1 and 2) 
since the distribution of sources used to obtain the desired thickness 
distribution inevitably produces a camber distribution of its own. This 
“camber due to thickness” must then be included in the distributed 
vortex problem if one is to obtain the desired overall camber. 

The results given in this paper, however, will be restricted to those for 
the indirect lifting problem: given the loading, and omitting thickness 
effects, what is the associated camber line of the blades, how does the flow 
and pressure field develop, and what is the performance of the rotor? The 
last of these questions has been substantially answered in reference 9, 
since most performance characteristics calculable in this theory do not 
depend on the details of the chordwise loading distribution but primarily 
on r(r). This is not to say that the chordwise loading distribution is not 
important; for example, it will affect boundary-layer behavior, stall 
margins, etc., of the rotor and is of great interest for these and other 
reasons. 


RELATION TO AXISYMMETRIC “THROUGH-FLOW” 
THEORY 

Let us consider first the mean pressure level (i.e., the azimuthal average 
pressure) developed by the rotor. We have generally, in the strictly 
linear approximation, 

p=p— (p_ OJ )= —p-aoU\/l+<T 2 V x - 

where 

v x > = v x cos <p+vt sin <p 

la / d<f> d<j>\ 

= uVi+^ 2 \dz + de) 

Then for any axial station downstream of the rotor x>c ax (max) , using 
(3) -(6) and averaging over 9 , we obtain 

<p(x,r) 

/d<t> d(f>\ 

= ~ WP - \dz + M/ 


( 12 ) 


(13) 
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Upstream of the rotor (x<0) the corresponding result is 


(14) 


<: P{x,r)) u -{V- K ) 

(iiB 


— tf 

4tt/3 2 4 


exp 




(:-r) 


Ro(i<okv)yok(^) (15) 


If we recall that all yok’s are zero when 7/cos <p is independent of t) (imply- 
ing “constant work” design, r = constant = T), we see from (15) that for 
a constant work rotor, there is no change in mean pressure (starting 
from upstream) until the rotor itself is reached, while from (14) we see, 
for the same design, that the entire mean static pressure rise is achieved 
within the rotor passage and the “far-field” value (compare eqs. (14) and 
(11)) is already attained at x = c ax <m ax) . 

On the other hand, if r(r) ^constant, then the y 0 *’s do not vanish and 
there is essentially exponential approach, away from the rotor, to the 
respective upstream and downstream values of the mean pressure level. 
This type of result is typical of axisymmetric theories of axial compressor 
flows (ref. 12), but the rate of approach is, of course, sensitive to the 
assumed area distribution of the flow annulus (red. 13). For the particular 
case we have considered ( ruj >'t = constant) , the “decay length” for 
approach to the upstream and downstream values is of order j3(r T — r H ) /x 
(seeeq. (7) and refs. 11 and 12). At large (subsonic) axial .Mach numbers, 
approach to the asymptotic states is very rapid, as expected (ref. 12). 

The result (1) or (14) indicates a constant radial static pressure 
profile far downstream of the rotor, at obvious variance with the need 
for a radial pressure gradient to balance the centripetal acceleration asso- 
ciated with the induced tangential velocities (see eq. (18)). This is 
simply a result of the linearization used throughout the present theory 
(for example, eq. (12)). However, comparisons between appropriate 
azimuthal averages of the results from the present three-dimensional (but 
linear) theory and higher-order (but two-dimensional) “actuator disc” 
results are expected to suggest means of identifying and including the 
more important nonlinear effects. We have already shown (ref. 9) that 
the second-order calculations, used in that paper to compute losses, are 
consistent with “radial equilibrium”. It should be possible to include 
certain nonlinear effects, such as centrifugal effects, consistently in a 
modified three-dimensional theory. Further work in this direction is 
underway. 

The results (14) and (15) can be understood in terms of the mean 
(azimuthally averaged) stream-surface deflections associated with non- 
uniform loading. (Note that the wake functions themselves, depending 
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only on ( and rj, yield zero first-order pressure perturbations.) If T(r) 
increases from hub to tip, more work is done by the rotor at the tip than 
at the hub. This appears initially (just behind the rotor) as a higher 
pressure near the tip than at the hub. In fact, using (7) and evaluating 
(14) for x~c ax (mHX> with the assumption c„ x (max> / r T <£\, we find, imme- 
diately behind the rotor, 

(p(c^,r)) d - {p-*)=- (r(ij)+r) (i6) 

where we have used the first of equation (9) and also equation (4) . The 
higher (lower) pressure at the tip (hub) must relax to the constant value 
given by the linear theory far downstream of the rotor. Moreover, the 
azimuthally averaged flow is effectively subsonic if M< 1 (ref. 12). This 
means contraction of the outermost streamtubes and expansion of those 
near the hub. The stream surfaces must, therefore, be deflected outwards. 
In fact, we find 

<t t Y f ( r\ 

5(>?)<*>-con S t = UL ^J^ (17) 


provided T— r«r. These average radial stream-surface deflections are 
zero, as required, at q = h and q — 1 and are positive if F increases radially 
and negative if I' decreases radially. There is no (average) streamline 
deflection if r = constant = T. 

Associated with (14), (15), and (17) are certain azimuthally averaged 
tangential, axial, and radial velocity profiles. Downstream of the rotor 
(see also ref. 9) , these are given by 


(Vt) d = 


BT( V ) 

2irr 


(18) 


{Vx) d 


wB 

2irU/3 '• 
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r-/3 2 r(r?) 


i =0 r / ,\ 


(Vr) d = 


4ir/3C/ 


„ (max) _ . . _ 

-uB ^ r c « , / co A 


yok(li) Ro(KQky) j+ U (19) 

yok(£) Ro 1 (kum) ( 20 ) 


The actual velocity field derivable from (3)-(6), of course, is curl-free 
except at the B helical sheets (wakes) of concentrated shed vorticity, 
this vorticity having strength proportional to dY/d-q and being oriented 
so as to lie in the helical sheets representing the wakes. 
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By contrast, the vorticity of the (mean) downstream velocity field 
defined by (18)-(20) is 

S2<">= curl °) (21) 

Thus, since Qe {v V^ v) = wr/ZJ, an important effect of defining an “equiv- 
alent axisymmetric” flow (through averaging the more detailed three- 
dimensional theory) is to replace the original concentrated vorticity by an 
equal total amount of vorticity which is, however, distributed uniformly 
over the flow annulus, yet still oriented along the zeroth-order streamlines. 
This distributed downstream vorticity is an important feature of axisym- 
metric through-flow theories (ref. 12); results such as (21), (19) and 
(14) help to establish connection between the latter theories and that 
described here. 

It should be clear from the preceding that the three-dimensional (but 
linearized) potential theory contains many important elements of the 
axisymmetric theories and has, in addition, the capacity to describe 
azimuthal variations superposed on those results. Some interesting 
examples of the latter are given in the following section. 


THREE-DIMENSIONAL PRESSURE FIELD (LIFTING 
PROBLEM) 

The azimuthally-averaged pressure fields (14) and (15) are the same, 
regardless of whether the compressor is transonic or not (ref. 12). But 
the azimuthal variations about these mean levels are vastly different, 
depending on (1) whether or not the tip relative Mach number is super- 
sonic and (2) whether or not the rotor is uniformly loaded (along the 
blade span) . 

The latter observation comes from consideration of the conditions 
under which transonic “acoustic resonance” (see the introduction) can 
occur for the strictly lifting problem. As mentioned earlier, when 
M 2 rT = M i (l J r<TT i )> 1, some of the quantities X „4 (eq. (7) ), as determined 
by the linear in viscid theory, can vanish. A glance at equations (5) and 
(6) shows that when this happens some of the acoustic eigenmodes can 
be amplified indefinitely unless the corresponding fe B *(?) are identically 
zero. This is the “resonance” to which we have referred. However, we have 
already noted that the fi„*(£) are, in fact, zero if the span wise loading is 
constant at each f. Thus a constant work rotor with similar chord wise 
loading profiles (i.e., 7 (*?,£)/cos ip= Tg(£) ) excites no transonic reso- 
nance. A more general way of saying this is that any (purely lifting) 
transonic rotor which sheds wakes of free vorticity can excite acoustic 
resonances of spinning modes and otherwise not. (Finite values of h n k 
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can be regarded as representing excitation of acoustic modes by the 
wakes.) Of course, the finite thickness of the blade of any real rotor will 
always induce transonic resonances (ref. 2). 

If transonic resonance is excited, one or several acoustic modes will be 
singled out and tend to dominate the pressure field, although the (linear) 
azimuthal average of the pressure fields will vanish. Naturally, the 
resonant modes do not attain infinite amplitudes; their finite amplitudes 
are easily predicted by including either viscous effects (as in refs. 1 and 
2) or certain nonlinear effects (as in ref. 10) or both. In the following, we 
present typical numerical results extracted from reference 10; included 
are pressure-field results for a typical rotor (.6 = 40, 6 = 0.8) operating 
both subsonically (M rI . = 0.9) and transonically (M rT = 1.054). In refer- 
ence 10, both a uniformly loaded rotor (r = T) and a rotor with linearly 
increasing T(jj), characterized by an increase of approximately 20 percent 
in T from hub to tip, were analyzed. However, for lack of space we present 
here only the results for a uniformly loaded rotor. The axial Mach number 
in the examples discussed is M = 0.5. The specific loading distributions 
used in the calculations were of the factorizable type y/cos < p= r(r/)g(£), 
where 

ff(£) ; V£(c»x-£) 

irc a x z 

Note that this chord wise loading distribution is symmetric about the 
midchord, nonsingular, and satisfies the Kutta condition and the nor- 
malization (11). 

Our first example (fig. 3) shows the “near field” subsonic upstream 
pressure fluctuation, at x/c KX = —0.1, over an azimuthal period AO = 2ir/ 6 
(corresponding to a single blade passage). Five radial stations are indi- 
cated. The tip Mach number being 0.9, we have <jt= 1.46. The pressure 
side of the airfoil corresponds to the smaller values of B0/2ir. For the 
solidity noted in figure 3, the projection of the helix of advance forward 
from x = 0 (for which BO/ 2v=\ is the blade leading edge) to x /c ax = —0.1 
corresponds to 60/2ir = O.35. Note that the high- and low-pressure regions 
remain relatively well identified just upstream of the airfoil (typical of 
subsonic flow) . The fluctuating signal is superposed on a zero mean pres- 
sure level, in agreement with (15) for a uniformly loaded rotor. This 
pressure signal is almost entirely dominated here by the (n=l, A;=l) 
eigenmode, not because of any resonance effect, but because the exponen- 
tial decay of this mode (for 6 = 40) is significantly slower than that of 
the remaining modes. 

Let us contrast this behavior (fig. 4) with the same rotor of the same 
solidity, but at the transonic state of M, T — 1.054. In this case, <tt= 1.84, 
and the angle (Bd/2v) of projection of the leading edge is approximately 
0.3. The fluctuating pressure signal is significantly distorted, relative to 
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Figure 3. — Pressure variation just 
upstream of typical uniformly 
loaded rotor operating at subsonic 
tip Mach numbers. Angular range 
corresponds to one full blade 
passage. 


UNIFORM LOADING, D=0 

B - 40, h = .8, C ox /Lf 1.06, M=.5, M T = .9, <t t =I.46, X/C ox =-O.I 



UNIFORM LOADING, D=0 



- 0.1 


Figure 4. — Pressure variation just 
upstream of typical uniformly 
loaded rotor operating at super- 
sonic lip Mach numbers. Note dis- 
tortion due to modes above cut-off. 


the subsonic case, by the presence of a large number of modes which are 
now above cut-off and therefore propagate. (For example, the k= 1 modes 
for all n are above cut-off at this tip Mach number; the (7,2) mode is also 
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above cut-off, etc.) The shift toward smaller angles of the expansion region 
just ahead of the blade is consistent with the expected formation (for the 
specified loading) of an expansion fan emanating from the leading edge 
of the blade tip (rj= 1), but the flow field is still predominantly subsonic 
in character over the whole annulus. Hence, no really clear-cut develop- 


ment oi a quasi-two-uiiiiensional type flow hold can be idem 


UUCU 111. 


region. 

A clearer distinction between hub and tip section contributions to 
the pressure pattern seems to develop just behind the rotor (figure 5) 
indicating some semblance of quasi-two-dimensional behavior. Here, the 
pressure signal oscillates about the uniform downstream level given by 
(1) or (14) at each ?j (because of the uniform loading) but a fairly strong 
local compression followed by an equally strong expansion fan seems to be 
emanating from the aft portions of the blade tips. 

The results available in reference 10 indicate similar behavior for the 
nonuniformly loaded rotor, except that the near-field upstream and down- 
stream pressure signals oscillate about different mean levels for each rj, 
in accord with (15) and (14). One other major difference occurs for the 
case of nonuniform r, however; namely, the expected transonic resonance 
appears (see earlier discussion). For Mtt= 1.054 the (1,1) mode is the one 
nearest resonance, and the relative fluctuations in the near-field pressure 
that it produces are noticeably larger than for the constant-r case. The 
radial structure of this resonant mode also makes it more difficult to 
separate hub and tip behavior in the strip theory sense. 

Far-field pressure signals (for example, at x/c ax =— 5.0, x/c a * = +6.0) 
were also computed in the work reported in reference 10. As expected, 
essentially no pressure fluctuation about the azimuthal mean is observed 


UNIFORM LOADING. D= O 

8 = 40, h= 8. C„/L t = I.06, M = .5, M T =I.054, o- T =l.04. X/C„,= I.I0 



Figure 5. — Pressure variation just 
downstream of the same rotor as in 
figure 4- Note development of a 
semblance of a quasi-two-dimen- 
sional pattern. 
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at subsonic tip Mach numbers, while substantial fluctuations, com- 
parable to the near-field values, are found to be present as soon as M TT 
exceeds unity. 

CONCLUDING REMARKS 

In this paper, we have discussed some of the aspects of the three- 
dimensional aerodynamic theory of an axial compressor rotor. Our main 
emphasis has been on the newly developed lifting-surface (distributed- 
vortex) theory which forms the complement to the previously developed 
(distributed-source) theory for blade thickness effects. 

In the third section, we showed that the present theory contains many 
aspects of the axisymmetric through-flow or actuator-disc theories, while 
in the fourth section we gave examples of the type of additional informa- 
tion that the three-dimensional theory offers. 

It may be recalled that in reference 1 the relationship between three- 
dimensional and quasi-two-dimensional cascade theory was emphasized. 
Not surprisingly, in the lifting problem just as in the thickness problem, 
cascade theory, and corrections to it, can be derived readily from the 
three-dimensional theory. The techniques involved are similar to those in 
reference 1 ; details and some examples are available in reference 10. 

We find, in short, for the purely lifting case, that quasi- two-dimensional 
cascade theory is an excellent approximation, at least within and near the 
blade passage, for uniformly loaded rotors, even in the transonic regime. 
However, for a rotor with only moderate variations in spanwise loading, 
we find that the wake-induced velocities have a surprisingly large effect on 
the effective incidence of the blades. Much work remains to be done, 
therefore, before the complete relationship between cascade theories 
(and/or data) and the aerodynamics of three-dimensional compressor 
rotors will be fully understood. Our expectation is that the present theory, 
despite the limitations imposed by our assumption of small disturbances, 
neglect of viscosity, etc., will be a useful key in relating all the various 
approximate theories to the full three-dimensional problem. It has pieces 
of each part of the puzzle. 
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DISCUSSION 


J. C. YRANA (McGill University) : Do you find that, as blade number 
is increased, the 3-dimensional solutions converge toward the infinitely 
bladed 2-dimensional approximations, regardless of circulation dis- 
tribution? Going now in the opposite direction (as I have been mostly 
involved with low aspect ratios) would you agree that below a certain 
number of blades (10 to 12) it becomes impossible to design them for a 
prescribed variation of circulation (say 20 percent variation in circulation 
from mean)? 

McCUNE (author) : With regard to the first question, the 3-dimen- 
sional solutions converge, provided you are not in the transonic regime, 
to the 2-dimensional cascade solutions in the limit of large blade number 
and hub tip ratio approaching unity. This we showed for the thickness 
case in 1956 (refs. 1 and 2) and for the lifting case in reference 9 of the 
paper. On the other hand, for the “infinitely bladed” 2-dimensional 
approximation, by which I take it you mean the axisymmetric through- 
flow theory emphasized in the present paper, the 3-dimensional solutions 
will converge to that 2-dimen.sional approximation as B — > oo , regardless of 
regime. These are, however, formal mathematical results; I believe that 
in almost any practical case 3-dimensional effects play a role. 

Going in the other direction, the answer would have to depend on what 
you mean by “impossible”. With present theory, insofar as it does not 
include induced velocity effects, I think it would be impossible to correctly 
design the blades, even for blade numbers as high as 40 or 80 (let alone 10 
or 12), with 20 percent variation in circulation from mean. This is espe- 
cially true in the transonic regime. While the theory presented in our 
paper is limited by the linearizing assumptions, it is intended as a step in 
the direction of making such design possible. We intend to pursue it. An 
important step is to include thickness effects (ref. 1) since they induce 
camber of themselves. 



